Probing Quantum Frustrated Systems via Factorization of the Ground State 
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The existence of definite orders with simple periodicity in frustrated quantum systems is related rigorously 
to the occurrence of fully factorized ground states below a threshold value of the frustration. Ground state 
separability thus provides a natural measure of frustration; Strongly frustrated systems are those that cannot 
accommodate for unentangled and mean field solutions. The exact form of the factorized ground states and the 
critical frustration are determined for various classes of non-exactly solvable spin models with different spatial 
ranges of the interactions. For weakly frustrated systems, the existence of disentangling transitions determines 
the range of applicability of mean field descriptions in biological and physical problems such as stochastic gene 
expression and the stability of long-period modulated structures. 
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Introduction.- Interest in frustrated quantum systems is due 
to the fact that they exhibit a large ground state degeneracy, 
or quasi-degeneracy, associated to complex structures of the 
quantum phase diagrams |l]]. Moreover, frustrated spin mod- 
els arise naturally in a variety of physical situations, e.g., in the 
study of high-Tc superconductivity iSS, long -period modu- 
lated structures of condensed matter 141], and biological sys- 
tems with stochastic components ijstl ■ Unfortunately, rigorous 
results are available so far only for ultra-simplified models jit], 
while numerical simulations are challenging, because quan- 
tum Monte Carlo methods are not practical for frustrated spin 
and fermion models 01 . and the density matrix renormaliza- 
tion group is difficult to apply to systems with dimensionality 
larger than one and/or periodic boundary conditions ^ . Fur- 
thermore, criteria to quantify frustration and the classification 
of "weakly" and "strongly" frustrated systems have been so 
far phenomenological and/or numerical, with little or no rig- 
orous physical and mathematical insight l^l. 

Recently, a method based on quantum information tech- 
niques 121, has allowed us to establish rigorously that large 
classes of (generally non exactly solvable) frustration-free 
quantum spin models admit totally disentangled ground states 
(GS) at finite values of the interaction strengths and of the ex- 
ternal fields ifioll . These factorized GS coincide exactly with 
the mean field solutions and identify well defined magnetic 
orders endowed with simple spatial periodicities. In this Let- 
ter we build on that formalism to investigate the physics of 
quantum frustrated systems by showing that GS separability 
is a measure of quantum frustration that discriminates quan- 
titatively regimes of weak and strong frustration. Indeed, we 
prove the existence of a critical frustration (frustration thresh- 
old) below which fully factorized GS are allowed and corre- 
spond to a definite magnetic orders with simple periodicity 
for quantum spin models belonging to different universality 
classes and with different types and ranges of interactions. 
We then prove rigorously the existence of disentangling tran- 
sitions in the GS in the regime of weak frustration and identify 



the different magnetic orders and their quantum phase bound- 
aries. That, indeed, entanglement and separability can be used 
to qualify and quantify frustration (and viceversa) can be in- 
tuitively understood by observing that the presence of frus- 
tration tends to enhance correlations among the constituents 
and thus to depress the possibility for the occurrence of sep- 
arable (uncorrelated) states. We also discuss some predictive 
consequences of these results on the study of complex sys- 
tems in physics and biology. Although in what follows, for 
ease of presentation, the analysis is carried out in detail for 
one-dimensional systems, it is straightforward to extend it and 
apply it in general to systems of arbitrary dimensionality. 

Consider, without loss of generality, spin-1/2 systems for 
which frustration arises from the simultaneous presence of 
competing antiferromagnetic exchange interactions of differ- 
ent spatial range. The anisotropy > {a — x, y, z) in 
the spin-spin coupling at sites i and j of the lattice with dis- 
tance r — |i — j| is taken independent of r, and all couplings 
are rescaled by a common, distance-dependent factor fr > 0. 
The general model Hamiltonian then reads 

H =y~^/r {JxSf Sfj^f^Jy Sf Sf_^_^J 2 S'f S^_^_^) — h'^^ S^, ( 1 ) 

z,r<rniax i 

where Sf are the spin-1/2 operators at site i; h is the ex- 
ternal magnetic field; and rmax > 1 is the interaction range, 
i.e. the maximum distance between two spins with nonvanish- 
ing coupling. Without loss of generality we assume > Jy. 
This general XYZ Hamiltonian with interactions of arbitrary 
range includes as sub-cases very many different models span- 
ning several classes of universality, such as the short- and 
long-range Ising, Heisenberg, XY, XX, and XXZ mod- 
els. We first recall briefly the basic findings on GS fac- 
torization in frustration-free spin models jlOtl . The quan- 
tity controlling GS factorization is the entanglement excita- 
tion energy (EXE) AE |@]. At every site k it is defined as 
AE = iniinu,}{'^\UkHUk\G) - {G\H\G). Here |G) is 
the GS of the system and Uk is any local rotation acting on 
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the spin at site k, i.e. a single-spin unitary operation 1191]: 
Uk = ®i^k 1* ® "^Ok, where 1^ is the identity operator on 
all the spins but the one at site k, and Ok is a generic Hermi- 
tian, unitary, and traceless operator For any translation- 
ally invariant and frustration-free Hamiltonian H such that 
[H, Uk] 7^ VC/fc, the vanishing of the EXE is a necessary 
and sufficient condition for GS factorization In fact, the 
minimization defining the EXE identifies an extremal opera- 
tion C/fe at each site of the lattice and, therefore, a global op- 
erator U = ^i^Uk that admits as its own eigenstate the fully 
factorized (separable) state \Gf) 



is: 



-J, 



y -r COS 9Jx + COS Lpk COS ipk+r sin OJz = . (4) 



\Gp) = n [cos((?fc/2)| U) + e"""- sin(0fc/2)| 1^)] , (2) ''^ " 2^^^^^ \ /^XJ;; 



Because Eq. (HI must be satisfied both for the cases in which 
fk = Vk+r and when 7^ "y^fc+r, it must be either sin6' = 
or Jz = 0. The first case {sva.6 — 0) is trivial, as it im- 
plies saturation rather than proper factorization. The second 
possibility ( = 0) is associated with proper nontrivial fac- 
torization, characterized by 6* 7^ 0. Using Eq. (|4| and Eq. Q 
one determines exactly the factorizing field: 

r /< 1/2 



/ > 1/2 • ^ ^ 



where 9k and ipk define the direction of tJk in spin space. 
\Gf) is the exact GS if and only if the EXE vanishes. Proceed- 
ing to investigate frustrated spin models by this method, we 
consider first the simplest short-range interactions, i.e. anti- 
ferromagnetic interactions that extend up to nearest-neighbor 
(nn) and next-nearest-neighbor (nnn) spins (rmax = 2). 

Short-range models of frustrated antiferromagnets.- For 

''max = 2, with ordering /i > /2 > 0, in the Hamiltonian 
Eq. dTji we can set, without loss of generality, /i = 1, so that 
the parameter / = /2//1 S [0,1] quantifies the degree of 
frustration: For / = the system is frustration-free, while for 
/ = 1 the model is fully frustrated. To verify the existence 
of a factorized GS we impose minimization of the energy and 
the vanishing of the EXE to determine the explicit expressions 
of Ok and (fk- One finds that as long as /< 1/2 a factorized 
GS exists and is associated to the single-step antiferromag- 
netic (SA) order along the x axis. This behavior is mirrored 
in the fact that fk ~ kn, Vfc. Viceversa, as soon as / > 1/2, 
the candidate factorized GS is associated to a dimerized anti- 
ferromagnetic order (DA), corresponding to alternating local 
phases: ip2k = kir, ip2k+i = f-ik (See Fig.!!). The angle 9k 
is site-independent: Ok = Vfc, and is the solution of 
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(3) 



where hp stands for the factorizing field, i.e. the value (at this 
stage, yet to be determined) of the external field at which the 
GS becomes fully separable. The quantities JTq are the net 
interactions that express the coupling of the entire system to 
a given spin, due to the presence of the external field. The 
net interaction along the z-axis is independent of the mag- 
netic order: Jz = 2J2(1 + /), while for the ones along a; 
and y one has Jx.y = ^2(1 — f)Jx.y in the presence of 
SA order (/ < 1/2), and Jx^y = —2fJx,y in the case of 
DA order (/ > 1/2). To prove that the state in Eq. Q 
is an eigenstate of the Hamiltonian we decompose the lat- 
ter at h = hp into a sum of terms involving only pairs 
of nn and nnn: Hk,k+r = fr{JxS^S^^.,, + JyS^Sl^.,, + 
Jz S^Sl^^ ) - /ij (S'f _|_^ ) , where, consistently with Eq. Q, 
^/ = frCOsO\,Jz—cos{Lpk)cos{tfk+r)Jx]/'^- Thus the Condi- 
tion for \ Gf) to be an eigenstate of every pair interaction term 



A sufficient condition for \Gf) to be the GS is that its projec- 
tion over every pair of spins be the GS of the corresponding 
pair Hamiltonian iflolfTin . This condition is never satisfied in 
the presence of frustration, whose effects cannot be captured 
by quantities involving only pairs of spins. The method must 
be generalized to include minimal finite subsets of spins en- 
compassing frustration. In the case of rmax = 2, the minimal 
subset is any block of three contiguous spins, tagged fc — 1, 
fc, and fc + 1. The corresponding triplet Hamiltonian term 
Hk = \Hk-i.k + \HkM+i + Hk-i.k+i includes all the dif- 
ferent types of irreducible interactions appearing in the model. 
Exactly at h = hp we have that H = J^k ^k- Moreover, the 
projection of \Gp) over the Hilbert space of the three spins 
fc — 1, fc, and fc + 1 is an eigenstate of Hk- Therefore, if one 
can show that the projection of |Gi?) is the GS of every three- 
body term Hk, factorization of the total GS is proven. The 
analysis yields that: (i) if Jy = 0, the factorized state Eq. ^ 
is the GS of the systems at h = hp for all values of the frus- 
tration / e [0, 1]; (ii) if Jy 7^ 0, the GS is factorized when / 
lies below a critical value fc'. 
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yJxJy + Jy 
Jx Jy 



(6) 



To assess whether for f > fc there may be still a region com- 
patible with GS factorization we consider a partition of the 
Hamiltonian into blocks of more than three spins. We define 

the sequence of operators (i?^"'' = 'J2j=-n-^>'+'y ) which, 
for any integer n, admit the {2n + l)-spin projection of |Gi?) 
as their eigenstate, and whose lowest eigenvalue, in the limit 
of large ji, coincides with the GS energy of the total Hamil- 
tonian H. For every n, the eigenvalue of i/^") associated 
to the factorized eigenstate is e(n) = (2n — l)£p, where 
Ep is the energy density per site at h — hp. Denoting by 
/i(n) the minimum eigenvalue of i/*^"', we have that only 
if there exists an integer n such that A(n) = /i(n) — e(n) 
vanishes for any n > h, then the factorized state is associ- 
ated to the lowest eigenvalues of i/^"', and hence it is the 
GS of the total Hamiltonian H. By studying A(7i) as a func- 
tion of n one can determine exactly, albeit numerically, the 
actual boundaries separating the occurrence and the absence 
of GS factorization, as reported in Fig. [T] The exact thresh- 
old value ft lies just slightly above the analytical lower bound 
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FIG. 1: (Color online) Analytical lower bound fc (solid blue line) 
and exact numerical value of the frustration compatibility threshold 
ft (dashed black line) as functions of the ratio Jy/ Jx- GS factoriza- 
tion occurs if and only if / < ft. For / below the horizontal dotted 
line / = 1/2 the magnetic order is single-step antiferromagnetic 
(SA), while for / above it, it is dimerized antiferromagnetic (DA). 



Therefore, no factorized GS supports DA order, except at Jy 



0. 



fc, Eq. (|6]l. Summarizing, we have shown that for / < ft, 
short-range frustrated models admit an exact, separable GS 
of the form Eq. (|2]i when the external magnetic field h takes 
the value hp defined by Eq. (|5]l. The associated magnetic 
order is SA, while every factorized state associated to a DA 
order is always an excited energy eigenstate; it becomes a 
GS only exactly at Jy = 0. Because the existence of fac- 
torized energy eigenstates is always associated to a violation 
of the parity symmetry of the magnetization along the direc- 
tion of the external field [Till the proof of the existence 
of factorized GSs, even if it yields no direct information on 
the location of the quantum critical points he, warrants the 
existence of quantum phase transitions to an ordered phase, in 
frustrated models, as the external field h decreases and crosses 
he (in the case that we have illustrated, it is a transition to 
a SA order along x). Moreover, since the factorizing field 
hp necessarily lies in the ordered region, one can, at least, 
conclude that the factorizing field anticipates the critical one 
from below; hp < he, a behavior already evidenced in some 
frustration-free models lll2ll . Exactly at f = ft the system 
undergoes a level crossing, and hence a first order phase tran- 
sition from the twofold degenerate factorized GS Eq. (|2]i to 
a twofold degenerate entangled GS state with complex long 
range incompatible with factorization points because the lat- 
ter are necessarily associated to ferromagnetic or antiferro- 
magnetic orders along a fixed axis ifioll . This phenomenon 
identifies a frustration-driven entangling-disentangling transi- 
tion of the GS at h = hp as f crosses the critical threshold ft 
dividing the regimes of weak and strong frustration. 

A phenomenological measure of frustration in antiferro- 
magnetic models is provided by Tcw/Tn i-e. the ratio of 
the Curie-Weiss temperature to the Neel temperature of bulk 
three-dimensional ordering lll3ll . This definition cannot be ap- 
plied to systems with vanishing Tn (like, e.g. 1-D and 2-D 
models) and it cannot distinguish between different contribu- 
tions, classical and quantum, to frustration. Notwithstanding 



these limitations, our ground-state analysis suggests that this 
phenomenological measure does capture some aspects of frus- 
tration, as follows. The existence of a factorized GS implies 
the existence of a transition to an ordered phase (occurring at 
a critical field higher than the factorizing field he > hf). In 
a bulk 3-D system this should coiTespond to a finite value of 
the ratio Tcw/Tn, since the order would necessarily freeze 
at a finite value of Tn > 0. On the other hand, in strongly 
frustrated models, strong correlations between quantum fluc- 
tuations persist in the presence of an applied field favoring 
GS factorization, and the existence of ordered phases tends 
to be suppressed. In the corresponding bulk 3-D systems one 
should then find a higher, asymptotically diverging value of 
the ratio Tcw/Tn, as the temperature at which the order 
freezes approaches the absolute zero. 

Models with interactions of arbitrary finite range.- For 

models Eq. ^ with finite J'max > 2 the triplet Hamiltonians 
Hk are generalized to subsets of r,„ax + 1 spins, with con- 
straint J2k = H ath = hp. The space of the Hamiltonian 
parameters is still divided in a region of low frustration com- 
patible with GS factorization, and one of high frustration for 
which GS factorization is forbidden, as shown in Fig. |2]for 
models with maximum range of interaction r„iax ~ 4 that in- 
clude, for instance, the ,/i — J2 — -h ~ Ji and Ji ~ J2 — J3 
models. Two general trends are observed: For Jy ^ 0, the 
factorized GS has SA order along the x-axis and, as shown in 
Fig. [T] and Fig. |2l the region of low frustration allowing GS 
factorization decreases as the anisotropy Jy/ Jx increases. 

Models with interactions of infinite range.- If in Eq. ([T]) we 

let fr when the maximum interaction range rmax ^ 
00, the question of the existence of factorized energy eigen- 
states can be analyzed by neglecting all the interactions be- 
tween spins at distances greater then some cutoff value r' , 
solve the associated constraints, and then let r' 00. To 



this aim, for each r' we consider the operator 



(r') 



2 ^7=-r' ^7'=-r'v^ '^7-7W 2r'+i- 17-7' | -^7,7' that ex- 
presses the sum of all the pair interaction terms between the 
r' spins closest to k, and the associated quantity A(r') = 

l^i'r') + \{Jx + Jy) TJk=i{-'^f fk, where /i(r') is the lowest 
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FIG. 2: (Color online) Threshold value of the frustration, ft, below 
which GS factorization occurs, for frustrated antiferromagnets with 
'"max = 4, as a function of Jy/Jx. Solid black line: /2 ~ f', fs = 
f; /4 = f\ Dashed red line: /a = f; fi = f/2; fi = //3. 
Dotted blue line: /z = f; fs = //2; fi = //4. 
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fr') 

eigenvalue of iJ j, . We have then analyzed different decay 
laws for fr, i.e. fast: fr = intermediate: fr = and 
slow: fr ~ l/'vA'- '^he first case it is always A(r') = 
and consequently the system admits a factorized GS exactly 
at hp = (7r^/12) y^JjTjy. In the second case (/,. = 1/r), 
according to the numerical evidence, A(r') vanishes in the 
limit of arbitrarily large r' and GS factorization appears to 
occur at hp = \n{2) J^Jy Finally, in the case of slow de- 
cay {fr = l/\/r), one has that A(r') ^ OVr' and therefore 
no factorized GS can exist. Therefore, fully connected mod- 
els characterized by a rapidly decaying fr, and hence by low 
frustration, allow for GS factorization and the associated SA 
or DA orders. Viceversa, models with slowly decaying fr, 
corresponding to strong frustration, do not admit factorized 
GS and simple mean-field and classical-like descriptions. 

Frustrated quantum models of complex condensed matter 
and biological systems.- The ANNNI (Axial Next-Nearest- 
Neighbor Ising) model, a particular case of the general class of 
models that we consider in the present work, provides a pos- 
sible effective description of systems with long-period modu- 
lated structures, such as, e.g., polytypism, anti-phase bound- 
aries in binary alloys, and helical phases in rare earths com- 
pounds I'Tl. It is thought that quantum frustration effects may 
be the mechanism responsible for the observed stability of 
these structures, and for this reason the quantum version of 
the ANNNI model is being intensively studied 101. It is then 
important to establish whether stable modulated structures are 
indeed predicted at all by the quantum ANNNI model and in 
what physical regimes. When applied to the quantum ANNNI 
model {Jy = = 0), our analysis proves that the mean- 
field description is applicable for all values of the frustration 
and that the value f ^ 1/2 discriminates between two types 
of stable structures, a simple unmodulated ferromagnetic or- 
der associated to a fully factorized GS for / < 1/2 and an 
anti -phase modulated GS with DA order for / > 1/2. 

Models of frustrated quantum spin networks have also been 
advocated as effective descriptions of gene expression and 
complex genomic patterns fsl. Here, again, the problem arises 
of the range of applicability of simple mean-field descriptions 
corresponding to simple magnetic orders. Indeed, much as for 
neural networks, the landscape of stable attractors in gene net- 
works depends, classically, on the degree of frustration. As- 
suming a description based on frustrated classical models with 
long-range interactions leads to the qualitative prediction of a 
small number of stable attractors in the presence of a "suf- 
ficiently" weak frustration. The question is then whether this 
prediction is stable against the effects of quantum fluctuations. 



Our analysis shows that the mean-field picture is qualitatively 
correct and makes it quantitative by determining the boundary 
between the weak and the strong frustration regime, in which 
the mean-field predictions fail. For models with long-range 
interactions, as we have shown above (see Fig. Fig. |2]i, the 
region of low frustration consistent with a mean-field descrip- 
tion is determined by the anisotropy ratio Jy / and decreases 
as the latter is increased. 

Conclusions and outlook.- We have introduced a rigorous 
criterion for discriminating between weakly and strongly frus- 
trated quantum systems in terms of GS factorizability. We 
have determined the threshold that separates the regions of 
weak and strong frustration, and we have singled out the ex- 
act forms of the factorized GS, the associated quantum phases, 
and the corresponding magnetic orders in the region of low 
frustration. These criteria should be experimentally testable 
in two- and three-body correlation experiments. 
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